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Intrinsic nonlinear response of surface plasmon polaritons 
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We offer a model to describe the intrinsic nonlinear response of surface plasmon polaritons (SPPs). 
Relation of the complex nonlinear coefficient of SPPs to the third-order nonlinear susceptibility of 
the metal is provided. As reported in a recent study, gold is highly lossy and simultaneously highly 
nonlinear due to interband absorption and interband thermo-modulation at a wavelength shorter 
than 700 nm. The effect of the high loss of the metal on the SPP nonlinear propagation is taken 
into account in our model. With the model we show difference in sign of real and imaginary parts 
between the nonlinear propagation coefficient and the nonlinear susceptibility of component material 
for the first time to our knowledge. Our model could have practical importance in studying plasmonic 
devices utilizing the nonlinear phase modulation and the nonlinear absorption of SPPs. For example, 
it allows one to extract the complex nonlinear susceptibility of gold through a measurement of SPP 
nonlinear propagation at the visible range. 

PACS numbers: 42.65.Tg, 73.20.Mf, 42.65.Hw, 42.65.Wi 


I. INTRODUCTION 

Nonlinear optics and nonlinear fiber optics have led 
to realization of modern photonic functionalities such 
as optical frequency conversion, generation of ultrashort 
pulses, all-optical signal processing and ultrafast switch¬ 
ing 00- In the last decades, many efforts have been 
to the study of nonlinear optical processes assisted by 
plasmonic nanostructures 0 because of nanoscale field 
confinement, field enhancement at metal-dielectric inter¬ 
faces 01, a strong and ultrafast third-order nonlin¬ 
ear response of metal 1 and other novel phenomena 
not achievable in conventional structures, for example 
extreme-ultraviolet continuum generations 0 ■ 

For extended metal-dielectric interfaces, surface elec¬ 
tromagnetic waves can propagate at the metal-dielectric 
interfaces, which are called surface plasmon polartons 
(SPPs). SPP-mediated nonlinear processes have been 
extensively studied, for example harmonic generation Q, 
four-wave mixing , nonlinear absorption [10J, self-phase 
modulation 11 -14] and soliton 15-19] and all-optical 
modulation 20]. 

Many theoretical studies of the nonlinear response of 
SPPs in different structures of plasmonic waveguides 
have been performed 11-H9I. 121-25]. A common theoret¬ 
ical approach to describe the nonlinear response in op¬ 
tical waveguides and fibers is based on the scalar model 
[§] disregarding influence of a longitudinal component of 
electric field which is justified for large core waveguides. 
Full vectorial model [27], El is developed for describ¬ 
ing the nonlinear response in subwavelength waveguides. 
In El, [l| the authors offered approaches to describe a 
surface-induced enhancement of the nonlinear response 
in plasmonic planar waveguides El and rod waveguides 
E|, which are in good agreement with the full vecto¬ 
rial model El . In most of these studies EMI , a linear 
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metal and a nonlinear dielectric were assumed and an 
imaginary part of permittivity of the metal e" is assumed 
to be small and perturbative. However, in practice the 
metal shows a strong nonlinear response 0, and SPPs 
are intrinsically nonlinear even in the absence of a non¬ 
linear dielectric. 


The nonlinear response of SPP based on the pondero- 
motive nonlinearity of metals in the infrared spectral 
range [HJ and the second-order surface plasmon solitons 
based on__the non-local multipole nonlinearities in metal 
surface E| were shown. 

Marini et al. E| studied the effects of the thermo- 
modulational interband nonlinearity of gold on the prop¬ 
agation of SPPs guided on gold nanowires at the wave¬ 
lengths longer than 750 nm. They also assumed low loss 
of gold which is justified at the long wavelengths. How¬ 
ever, at shorter wavelengths gold has a high loss [31] and 
thus an effect of the high loss of gold is needed to be con¬ 
sidered for modeling SPP nonlinear propagation. More¬ 
over, the quite large third-order nonlinear susceptibility 
of gold at these shorter wavelengths, as recently reported 
in [20 El, gives practical importance of modeling the 
SPP nonlinear propagation in the spectral range 1251. 12 til . 
De Leon et al. |25| developed an approach to describe 
the intrinsic nonlinear response of SPPs in planar waveg¬ 
uides taking account of the significant loss of gold e". 
They assumed only that the SPP loss coefficient is not 
too large. It is worthy to note that the condition of small 
SPP loss coefficient is not equivalent to the condition of 
small and even in the case of large e" the SPP loss 
coefficient can be small if the power portion in metal is 
small or the the SPP loss is compensated with a gain 
media in dielectrics [30 . However, this approach can be 
applied in planar waveguides, thus a general approach is 
needed. 

In this paper we offer a model to describe the intrin¬ 
sic nonlinear response of SPPs in a spectral range where 
the plasmonic metal is highly lossy, which is a modified 
version of the full-vectorial model. Our model shows a 
prominent difference in phase angle (one can find a hint 
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of the difference in phase angle in Fig. [3] of j25j) and 
even in sign between the complex third-order nonlinear 
susceptibility of the metal and the complex effective non¬ 
linear coefficient of SPPs. In particular, a positive imag¬ 
inary part of the complex nonlinear susceptibility of the 
metal, corresponding to a positive nonlinear absorption, 
can result in a negative nonlinear absorption (saturated 
absorption) of SPPs. We note that our approach can 
be used for all typ es of waveguides including the planar 
waveguides of [25j. 


II. THEORETICAL MODEL 

In this paper we assume the continuous wave and 
all the fields can be expressed as F (r, t) = (1/2) • 

F (r) exp (—iijot) + c.c , where c.c. signifies the com¬ 
plex conjugate. Below, we consider the time-independent 
fields F (r). For a linear plasmonic waveguide, we can 
write the electric and magnetic field vectors 

e (r) = \JZ 0 /s 0 * 0 (z) exp ( ikz) e 0 (rl), (1) 

h(r) = (iuifj, o) -1 V x e = 

\/l/(Z 0 So) v f , 0 ( 2 ) exp (ikz) ho (rl). (2) 


Here £ 0 , p .0 and Zq = is the permittivity, 

permeability and wave impedance in vacuum, respec¬ 
tively. z and z are a coordinate and the unit vec¬ 
tor in the direction of propagation and f± is a posi¬ 
tion vector in the transverse plane, so is defined as 
so = (1/2) / Re ^eo x hg) ■ zda, where the integral is per¬ 
formed in the transverse plane. The function To (z) has 
the form To (z) = To (0) exp (—az/2) and /? = n + ia/2 
is the SPP propagation constant. We note that the To 
is normalized so that the |To | is equal to the power flow 

along the 2 direction S (z) = (1/2) f Re x h*) ■ zda = 
|4-o| 2 . 

Now, we add the small perturbation of permittivity 8s 
to the linear plasmonic waveguide. If we assume a trans¬ 
verse distribution of the perturbation 8s = 8s (rl) and a 
longitudinal invariance of the perturbation d(8s) /dz = 
0, the perturbed electric and magnetic fields can be ex¬ 
pressed as follows. 


E (r) = v%/s 0 T (z) exp ( ikz) 
x [(eo (rl) + Sib (rl))], 

H (r) = () 1 V x E = 


( 3 ) 


Vl/^oso^ (z) exp (ikz) (h 0 (rl) + 8h 0 (rl)) . (4) 
Let’s remind of the Lorentz reciprocity theorem fl7l . 


d_ 

dz . 


Ei (r) xH 2 (r)- -^2 (r) x Hi (r) 


J [£2 (r) - si (r)] Ei (r) • E 2 (f) da, (5) 

where (^Ei,Hi) and (^E 2 ,H 2 ) are solutions of the 

Maxwell equations corresponding to the relative per¬ 
mittivity £1 (r) and s 2 (r), respectively. Let’s choose 

[Ei, Hi) and (e 2 , H 2 ) to be the unperturbed backward 
propagating field (^e~,h~) and the perturbed forward 

propagating field (jd , h). In the first order of perturba¬ 
tion, Eq. (]5]) leads to the following equation describing 
the amplitude T ( z ). 


dT a 

— = T + ik 0 ■ 5n ef f T, 
dz 2 


8 tl off — 


J 8s ■ ^e"o 2 — 2eg-^ da 
2^o / ^eo x h 0 ) 


■ zda 


( 6 ) 


( 7 ) 


Here we considered e 0x = e 0x , e 0y = e 0y , e 0z = —e 0z , 
h 0x = —ho x , h 0y = —hoy, /ip, = hoz. 

We emphasize that in the derivation of Eq. © and 
0 from Eq. © ©, no additional assumption has been 
taken, except for the small perturbation and longitudi¬ 
nal invariance of d£(rl). In particular, the assumption 
of small imaginary part of metal permittivity have never 
been taken, contrary to in the full-vectorial fl3l . fl4l . l27l . l28j 
and scalar [2j models. In our model the backward prop¬ 
agating fields are used as the multiplied ones and conse¬ 
quently the product terms on the left-hand side of Eq. © 
experience no loss and the imaginary part of metal per¬ 
mittivity on the right-hand side also disappears, while 
in the full-vectorial and scalar models conjugated fields 
were multiplied. We also note that in |l7j a formalism 
utilizing backward propagating fields was provided. 

To justify the model (7), we utilize this model in the 
single gold-air interface (e = £ go id • 9 (—a;) + d£ air ■ 9 (x)) 
for 8s = 8s 1 • 9 (—x) + 8s 2 • 6 (x) , where 9 ( x ) is the 
Heaviside step function, and compare with an analytical 
result, the full-vectorial model and the scalar model. For 
calculations we used the measured wavelength-dependant 
permittivity of gold jib] . If we use the formula of effective 
refractive index n eff = sjSiS 2 / (£1 + s 2 ) in the single 
interface |29j, we can obtain the analytical formula of 
Sn eif as following 


8n ei f — 


(£1 + d£i) (s 2 + ^£ 2 ) 


(£1 + fei) + (£2 + 8 s 2 ) Y £1 T £2 


£i£2 


( 8 ) 


By the full-vectorial model [HI, |27| and the scalar 
model [2j, 8n e ff can be expressed 


• zda = 


8tlq ff — 


f 8s ■ |e"o | 2 da 


2Z 0 J Re ^eo x h^) ■ zda 


(9) 
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FIG. 1. (Color online) Response of effective refractive index in the single gold-air interface to the perturbation of permittivity 
in the metal (a), (b) and the dielectric (c), (d), respectively. Real part (a), (c) and imaginary part (b), (d) of the response. 
The black line corresponds to our model ©, the blue circles to the analytical formula the red crosses to the full-vector 
model m and the green line to the scalar model 


^Tleff 


ko J 5e ■ |eo | 2 da 
2k J |eo | 2 da 


( 10 ) 


First we calculate for the case of fei = 0.01, 5e 2 = 0 
[Fig. © a) and (b)], where the perturbation is taken in 
the metal. As shown in Fig. [ 1 ] the result calculated by 
our model © coincides with that by the analytical for¬ 
mula © in all the considered range of wavelength. How¬ 
ever, it greatly differs from results of the full-vectorial 
m and scalar m models at wavelengths shorter than 
600 nm showing even difference in sign [see Fig. ©a)] be¬ 
cause the large imaginary part of permittivity of gold at 
these wavelengths influences a phase of the electric field 
in metal which results in difference in value of denomi¬ 
nator between Eq. © and Eq. © and (nnD- And our 
model shows an imaginary part of Sn eti produced by the 
real perturbation 5e which can not be described by the 
full-vector and scalar models [see Fig. ©b)]. 

For the case of 5e 1 = 0, §£2 = 0.01, where the pertur- 

I 


bation is taken in the dielectric, the result of our model is 
almost same as that of the full-vector model [see Fig.©c)] 
because even the large imaginary part of permittivity of 
gold doesn’t greatly influence the electric field in dielec¬ 
tric. 


III. INTRINSIC NONLINEAR RESPONSE OF 
SURFACE PLASMON POLARITONS 


To investigate an intrinsic nonlinear response of 
SPPs, we take the nonlinear perturbation Se = £ n i. 
If we take the self-nonlinear permittivity £ n i.seif = 

(3/4)y (3) (w;w,-w,u;) E 
lead to 


1.0, ,2f| , Eq. © and © 


dz 




( 11 ) 


(3/4) / X^ 3 ) |e"o 2 • 1 

(eS a 

- 2e 0z) 

| da 

Ml 

(eo x ho'j 

| • zda ■ f Re | 

(e 5 x h^j 

| • zda 


( 12 ) 


where 7 is the effective nonlinear coefficient and y® is 
the third-order nonlinear susceptibility. 


We apply this model [Eq. (fill) and (fT2l) ] to the single 
gold-air interface. It is assumed that the effective nonlin- 
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FIG. 2. (Color online) Effective nonlinear coefficient in the 
single gold-air interface according to the wavelength. The 
black line is by our model ©m the red crosses by the full- 
vector model m and the green line by the scalar model 00. 
(a) and (b) shows real and imaginary part of the effective 
nonlinear coefficient, respectively. 

earity results from the intrinsic self-nonlinearity of gold, 
while the dielectric is taken as a linear medium. The 
reported self-nonlinear susceptibility of gold x^ ©| is 
used. As expected, Fig. [ 2 ] shows that the effective non¬ 
linear coefficient by our model significantly deviates from 
those by other models at wavelengths shorter than 600 
nm. We can know from Eq. m and (USD that a phase 
angle of the complex effective nonlinear coefficient 7 by 
the full-vector and scalar model is same as one of the 
complex nonlinear susceptibility of gold x^- Therefore, 
it is clear from Fig. [2] that our model shows difference in 
phase angle between 7 and x*- 3 \ even sign-difference in 
real and imaginary parts of them. In particular, we can 
see a negative real part of 7 meaning a negative nonlinear 
phase-shift [see Fig. 0(a)] in spite of a positive real part 
of X < 3) at the wavelength of 535 nm and a negative imag¬ 
inary part of 7 meaning a negative nonlinear absorption 
[see Fig.©b)] in spite of a positive imaginary part of x^ 
at the wavelength of 502 nm. Our results are in moder¬ 
ate agreement with results of full-dimensional simulation 
(blue circles of Fig. [3]). 

As other examples of plasmonic waveguides, 64-nm- 
thick gold slab surrounded by air and gold rod waveguide 
with air cladding and the waveguide radius of 100 nm are 
calculated by our model. For calculations, the analyti¬ 
cal results for modes in metal slab ©| and metal rod 
[3] are substituted to Eq. m■ Fig.© a) and (b) shows 
the effective nonlinear coefficient in the gold slab (anti¬ 
symmetric mode) which is about ten times greater than 
that in the single gold-air interface. Fig. ©c) and (d) 




FIG. 3. (Color online) Power-dependence of phase shift for 
the wavelength of 535 nm and the propagation distance of 
1000 nm (a) and absorption for the wavelength of 502 nm 
the propagation distance of 500 nm (b) in the single gold-air 
interface. The black line is by our model, the red crosses by 
the full-vector model, the blue circles by a full-dimensional 
simulation. 

shows the effective nonlinear coefficient in the gold-rod 
(TM mode, m = 0) which is about billion times greater 
than in conventional fibers 0 and greater by several or¬ 
ders of magnitude than in nonlinear nanorod waveguides 
utilizing nonlinear dielectrics fl4j . The enhancement of 
nonlinearity in the examples of plasmonic waveguides is 
due to the field enhancement at metal-dielectric inter¬ 
faces HI and the strong third-order nonlinear response 
of metal, in particular the large interband nonlinear sus¬ 
ceptibility of gold ©| . 

It is expected that our model can be applied to ex¬ 
traction of a nonlinear susceptibility of metal from ex¬ 
perimental data. For example, in Ref. (25| , the authors 
developed a theoretical approach to describe the effective 
nonlinear coefficient in planar waveguides and extracted 
the third-order nonlinear susceptibility of gold x^ from 
a measured effective nonlinear coefficient in a gold slab 
at the wavelength of 796 nm which characterizes the 
nonlinear phase shift and nonlinear absorption experi¬ 
enced by the SPP ©I . For comparison, we obtained 7 = 
(1.16 + *1.00) x 10 ~ 7 V 2 /m 2 by our model in the single 
gold-air interface at A = 796 nm from the extracted value 
of X (3) = (4.67 + 13.03) x 10 ~ 19 V 2 /m 2 ©I , which is in 
ajjood agreement with 7 = (1.03 + *0.98) x 10 ~ 7 V 2 /m 2 

It is worthy to mention again that the accuracy of our 
model is limited by both the approximations used in our 
theory and in all the other models, the small perturba¬ 
tion and longitudinal invariance of Se. For validity of the 
first approximation, the power should be limited so that 
the nonlinear permittivity is small compared to the lin¬ 
ear permittivity. For a greater power, a nonperturbative 
nonlinear model [32 | should be employed. For validity of 
the second approximation, the condition that the effec¬ 
tive linear absorption a is small compared to the effective 
propagation constant /3 should be satisfied, which has 
been also used in other models including [Tj|, |3 [ 2 ^, |25| . 
This condition is satisfied if the SPP mode is appropri¬ 
ately differ from an exact resonance as in the cases above 
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FIG. 4. Effective nonlinear coefficient of antisymmetric mode in 60-nm-thick gold slab surrounded by air (a),(b) and TM 
mode (m = 0) in gold rod waveguide with air cladding and the waveguide radius of 100 nm (c),(d). The solid line is real part 
of the effective nonlinear coefficient and the dashed line is imaginary part of that. 


(Fig. OH]). 

It is predicted that if one uses a cross-nonlinearity with 
a pumping SPP at a longer wavelength instead of the 
self-nonlinearity, the total nonlinear phase shift and non¬ 
linear absorption can benefit from a longer propagation 
length of the pumping SPP at a longer wavelength and a 
large interband nonlinear susceptibility of metal [23| ex¬ 
perienced by signal SPP at a shorter wavelength. Based 
on the theoretical consideration [23| for the self-nonlinear 
susceptibility, the cross-nonlinear susceptibility of metal 
can be expressed as following. 

X (3) (w s ; w p , —Wp, <u s ) = -e 0 (u p £:" (cu p ) j T (w s ) • (13) 


If the pumping power greater than the signal power 
is introduced, the cross-nonlinearity is dominant in the 
total nonlinearity and Eq. ([6]) and © lead to 


(Mu 

dz 


= ——IE's 

2 


■ *27 c: 


4U m k 


(14) 


= ka 


(3/4) / \' (3) (w s ; w p , -Up, cu s ) 

eo K) 2 

eo 2 (cu s ) - 2 el z (u> s ) 

da 

ZoJ 

eo ((Us) x h 0 ((u s ) 

• zda ■ f Re 

^0 (^p) ^0 (^s) 

■ zda 


(15) 


where eo ( 04 ) and ho (cu k ) is the transverse distributions 
of electric and magnetic fields [see Eq. m and ©], re¬ 
spectively, at the frequency of 04 . The subscripts s and 
p denotes the signal and pumping SPPs, respectively. 
The maximum nonlinear phase shift and absorption is 
expressed as follows. 


Max (0ni) = 2Re (7 cross ) I’Up (0)1 /a, 


P’ 


Max In- 


I*- 


1®. 


=-4Im(7 cross ) |^ p (0)| /a v . (17) 


Finally, we would like to note that our model for the 
cross-nonlinear coefficient [Eq. (fl4l) and (fTSl) ] is fully ac¬ 
curate even at a SPP resonance for the signal, because 
a p is small and the longitudinal invariance of £ n i is fully 


(16) satisfied in spite of a large a s 
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IV. CONCLUSIONS 

We provided the relation of the complex SPP nonlin¬ 
ear coefficient to the third-order nonlinear susceptibility 
of the metal. The effect of the high loss of the metal on 
the nonlinear propagation of SPPs was fully taken into 
account. Backward propagating field instead of adjoint 
forward propagating one was utilized to derive the re¬ 
lation. We predicted that the high loss of gold due to 
interband transition can lead to difference in phase angle 
between the third-order nonlinear susceptibility of gold 
and the SPP nonlinear coefficient and even difference in 
sign of real and imaginary parts between them. In par¬ 
ticular, a negative imaginary part of the third-order non¬ 
linear susceptibility of gold, which normally corresponds 
to the negative nonlinear absorption or the saturated ab¬ 


sorption, can lead to a positive imaginary part of the SPP 
nonlinear coefficient, which means the positive nonlinear 
absorption or the reverse-saturated absorption. To our 
knowledge for the first time, the sign-difference between 
the nonlinear propagation coefficient and the nonlinear 
susceptibility of component material was reported. Al¬ 
though the high interband absorption of the metal can 
lead to a short SPP propagation length which wouldn’t 
be promising for nonlinear devices, the metal has the 
great thermo-modulational interband nonlinearity at the 
interband transition wavelengths [23l ]. Moreover, recently 
optical amplification of SPP by using gain media in com¬ 
ponent dielectrics has been demonstrated [33j and in¬ 
creasing the SPP propagation length could be a possible 
strategy to take advantage of the great interband nonlin¬ 
earity of the metal. 
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